Surface singularities in Eddington-inspired Born-Infeld gravity 
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Eddington-inspired Born-Infeld gravity was recently proposed as an alternative to general rela- 
tivity that offers a resolution of spacetime singularities. The theory differs from Einstein's gravity 
only inside matter due to nondynamical degrees of freedom, and it is compatible with all current 
observations. We show that the theory is reminiscent of Palatini f(7V) gravity and that it shares the 
same pathologies, such as curvature singularities at the surface of polytropic stars and unacceptable 
Newtonian limit. This casts serious doubts on its viability. 
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Recent years have witnessed a proliferation of alterna- 
tive theories of gravity, motivated by long-standing puz- 
zles in general relativity It is pertinent to develop 
theoretical benchmarks that will allow us to select the 
physically relevant candidates, and to subsequently use 
observational constraints in order to single out the few 
theories that are actually viable alternatives. 

One of the major riddles of general relativity is that it 
predicts the appearance of spacetime singularities origi- 
nating from regular initial data, e.g. in the gravitational 
collapse of massive stars and in the early universe. In or- 
der to resolve these singularities, an appealing proposal 
for a modified theory of gravity, the so-called Eddington- 
inspired Born-Infeld (EiBI) theory, was recently put for- 
ward in Ref. 0] and it has been subject to scrutiny in 
a number of works EiBI gravity is equivalent to 

general relativity in vacuum and does not propagate any 
degree of freedom other than a massless graviton. On the 
other hand, the theory introduces nonlinear couplings to 
the matter fields 0, 0] , which resolve at least some of the 
singularities appearing in Einstein's theory. 

The Big-Bang singularity in early cosmology is re- 
placed by a freezing or a bouncing behavior of the cosmo- 
logical scale factor, depending on the extra EiBI parame- 
ter [H . The gravitational collapse of noninteracting parti- 
cles does not lead to singular states in the nonrelativistic 
limit 0, Q . A tensor instability of the homogeneous and 
isotropic universe was found in Ref. Q and EiBI gravity 
has been also studied as an alternative to the inflation 
paradigm Possible constraints on the theory have 

been considered using solar models [4fl and cosmological 
observations dEH (see also Ref. However, a de- 

generacy between EiBI corrections and different matter 
configurations Q makes it difficult to put observational 
constraints without independent knowledge of the matter 
content of the theory. 

Previous literature on EiBI gravity mostly focused on 
phenomenological aspects of the theory; a more detailed 
study on its dynamics and on the structure of its field 
equations has not been performed yet. Here, we argue 
that the field equations of EiBI gravity have a peculiar 
differential structure which is similar to that of Palatini 



f(lZ) gravity [lfj and, as a result, they exhibit the same 
pathologies with the latter (see e.g. Refs. [lil - tlij and 
the review [HI) and with theories where matter is cou- 
pled to the Ricci scalar, which also have similar charac- 
teristics • Although these theories provide an appeal- 
ing (and in fact similar 21[) early-time cosmology, such 



pathologies cast serious doubts on their viability. 
EiBI gravity is described by the following action 0] 



S = 



d 4 x ( J | det (g a b + kTZ 



(ab)j 



(1 + kA)V? +S M [g a b^ 



M 



(1) 



where Sm [gab, ^ m] is the matter action, VP m generically 
denotes any matter field, lZ a b is the Ricci tensor built 
from the connection r^ b , g = | det(g a f,)|, A turns out to 
be the cosmological constant and k is the extra EiBI pa- 
rameter which has dimensions of length squared. Round 
(square) brackets denote (anti)-symmetrization. 

In the metric approach the field equations contain 
ghosts, which must be eliminated by adding extra terms 
to the action [22|, 23 1. Thus, EiBI gravity is naturally 
based on the Palatini formulation, i.e. the connection 
Y c ah is considered as an independent field. The origi- 
nal EiBI proposal is based on some crucial assumptions, 
namely: (i) the matter action is independent from r^ b ; 
(ii) the connection is symmetric, T c ab = Tg a ; (iii) only 
lZ( a b) appears in Eq. ([T]) and not lZ a b- This last assump- 
tion is often implicit. That is, it is common in the lit- 
erature to have lZ a b appearing in Eq. ([I}, even though 
any subsequent calculation is based on the implicit as- 
sumption 1Z\ ab ] — (or on the non-standard definition 
lZ a b = TZi a b))- This would be an extra constraint as, 



even for a symmetric connection, Tt[ a b] — ~®[bF a u and 
it does not vanish generically This extra constraint is 
not required provided that only the symmetric part of 
the Ricci is used in the action. See Ref. (24| for a similar 
discussion for generalized Palatini gravity. 

In principle, assumptions (i), (ii) and (iii) are not re- 
quired and, together with the action ([T]), they define 
a particular version of EiBI gravity. Relaxing (i)-(iii) 
would lead to a metric-affine version of the theory, sim- 
ilarly in spirit to the case of metric-affine f(TZ) theo- 
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ries [25|, but very different from the original theory. We 
will explore this possibility in a separate publication. 
Here we shall rely on the assumptions above, as in the 
original proposal [H . 

We start by expanding the action ([1]) at second order 
in k. This yields 



S 



1 



8irG 
+Sm [g a b, *m] 

where TZ = g ab TZ ab . 
such that 8ttG = 1. 
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For simplicity, we shall use units 
When k = 0, EiBI gravity reduces 
to the Palatini formulation of general relativity with a 
cosmological constant. As is well known, in this case the 
field equations impose that the connection is the Levi- 
Civita one and the theory reduces to Einstein's gravity. 
However, at first order in k, quadratic corrections in the 
curvature tensor built from the independent connection 
appear in the action @. The Palatini formulation guar- 
antees that, despite these extra terms, no higher deriva- 
tives of the metric field would appear in the field equa- 
tions. Note also that, when expanded order by order 
in k, the action ©takes the form of a specific Palatini 
f(R,ll ab ) theory [2l Hi . 

We now return to action ([1}. Independent variation 
with respect to the metric and the connection yields 



v c [Vot (o6) ]-v^^=o, 



(3) 
(4) 



where we have defined q ab 



g ab + nU {ab ) and V Q 
is the covariant derivative defined with T c ab , T ab = 
(g)~ 1 ' 2 SSM /Sg a b is the standard stress-energy tensor, 
whose indices are raised and lowered by g ab , whereas q ab 
is the inverse of q ab . After some manipulations, Eq. (U) 
takes the form 



1 



T ab = -^1 C (d a qbd + d b q ad - d d q ab ) . (5) 



On the other hand, using Eq. we obtain 
(1 + nA)g ab - nT ab 



q ab = 



y/g^dct ((1 + nA)g ab - nT ab ) 
which can be rewritten as 



(6) 



x [(1 + K A)g ab - nT ab 



g a b ■ 



(7) 



Equation ((5]) determines q ab algebraically in terms of g ab 
and T ab , whereas Eq. (|5|) determines T c ab as the Levi- 
Civita connection of q ab . Hence, one can use these equa- 
tions to eliminate T c ab from Eq. ([7]). Then, the latter is 
the equation that has to be used to determine g ab - 

It is straightforward to see that, after eliminating T c ab , 
Eq. ([7]) becomes a second-order partial differential equa- 
tion in g ab . However, it is equally straightforward to see 



that it also contains second derivatives of T ab . This is 
true in the full theory, but it becomes more explicit if we 
expand Eq. at first order in n. The expansion can be 
found easily by noting that 



ab 



ab 



where r ab = T ab — \g ab T + Ag ab . Using the expression 
above, we get a single equation for the metric g ab only: 



Rab— A<? & + T ab — -Tg ab + k 



Sab - Sg ab 



+ 2 [V a V b T - 2V c V (a r cb) + Ur ab } + 0(k 2 ) . (8) 

where S ab = T c a T cb — ^TT ab and we have used the fact 
that Tab is symmetric. Note that now R ab is built solely 
from the Levi-Civita connection of g ab and V Q is the co- 
variant derivative associated with g ab . Therefore, A in- 
side Tab does not contribute to the second line of Eq. ([8]) . 

General relativity corresponds to k = and for k ^ 
Eq. ([8J) contains second derivatives of T a b, i-e- at least 
third derivatives of the matter fields (unless we consider 
a fluid approximation of matter). This is in contrast 
to Einstein's theory, where usually only first derivatives 
of the matter fields appear on the right hand side of 
Einstein's equations. This different structure is also evi- 
dent in the Newtonian limit of the theory , which can 
be straightforwardly obtained from Eq. ([8]) but the re- 
sults holds for any value of k and does not hinge on 
the small k expansion. The modified Poisson equation 
|V 2 p, whose solution reads 



is V 2 $ = f 



$ = $ 



N 



4^ 



(9) 



where p is the matter density and $n is the standard 
Newtonian potential. Although the modified Newtonian 
regime has been studied in some detail and it provides 
interesting phenomenology H 0] , Eq. © shows that the 
gravitational potential $ is algebraically related to p. 
This demonstrates that gravity is noncommulative: un- 
like in Einstein's theory, the metric in EiBI gravity is not 
just an integral over the sources but it receives an alge- 
braic contribution from the matter fields and their deriva- 
tives. Any matter configuration which is discontinuous 
or just not smooth enough will produce discontinuities in 
the metric and singularities in the curvature invariants 
(which depend on the second derivatives of $), leading 
to unacceptable phenomenology. Clearly, this behaviour 
persist in the post-Newtonian limit (cf. Ref. [THJ] for a dis- 
cussion in Palatini f(lZ) gravity). These problems have 
been overlooked in the literature of EiBI theory 

As mentioned earlier, the same pathologies arise in 
Palatini f(TZ) gravity [15| and in theories with matter 
coupled to the Ricci scalar (joj . Indeed, the structure 
of Eq. © is the same as in Palatini f(lZ) theory (cf. 
Eq. (28) in Ref. [HI). A qualitative difference is that our 
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Eq. ((8j) contains derivatives of the full stress-energy ten- 
sor, whereas the field equations in Palatini f(1Z) gravity 
only contain derivatives of the trace T = g ab T ab . Thus, 
the problems we are discussing are already manifest when 
EiBI gravity is simply coupled to a Maxwell field, whose 
stress-energy tensor is traceless, T = 0. 

Let us now discuss some problems in construct- 
ing perfect-fluid equilibrium structures in EiBI gravity, 
which are related to its peculiar differential structure. 
We shall build on previous works which discuss similar 
pathologies in other theories 15, 2(|. Static and spheri- 
cally symmetric perfect-fluid stars were discussed in de- 
tail in Refs. [H [a]. Without loss of generality, a conve- 
nient ansatz for the metrics reads 

q ab dx a dx b = -p(r)dt 2 + h(r)dr 2 + rW, (10) 
g ab dx a dx b = ~F(r)dt 2 + B(r)dr 2 + A(r)r 2 dn 2 , (11) 

where we have used the gauge freedom to fix the func- 
tion in front of the spherical part of the auxiliary q a b 
metric. We consider perfect-fluid stars whose stress- 
energy tensor reads T a b = [p + P] u a Ub + g a bP, where 
u a = (1/VF, 0,0,0) and p(r) and P(r) denote the en- 
ergy density and the pressure respectively. 

Notice that the field equations © are simply alge- 
braic relations between q a b and g a b- Inserting the ansatze 
above into Eqs. one can solve for the coefficients of 
g ab in terms of q ab and T ab Q . Then, the dynamical equa- 
tions (|7|) do not explicitly contain second derivatives of 
the matter fields and can be solved for p, h and P, assum- 
ing an equation of state of the form P = P(p). Interest- 
ingly, in this formulation the field equations are equiv- 
alent to the standard Tolman — Oppenheimer — Volkoff 
equations for the metric q a b, with an effective equation 
of state 0- Since matter is covariantly coupled to the 
g a b metric, the standard conservation of the stress-energy 
tensor follows, V ' a T ab = 0. Finally, the interior solution 
is matched to the (unique) exterior Schwarzschild met- 
ric through appropriate junction conditions at the stellar 
surface [3, H|. However, in this formulation the physical 
g a b metric is not a dynamical quantity and a successful 
numerical integration does not necessarily mean that ge- 
ometric invariants, which involve derivatives of g a b, are 
regular. Indeed, we show here that the Ricci curvature of 
the physical g ab metric is divergent at the surface. This 
singularity has been missed in previous literature on EiBI 
gravity. 

For simplicity, we consider the asymptotically flat case, 
A = 0. Using Eqs. (fT0|) - ([TT|) . the Ricci curvature reads 

R g = [2r 2 A 2 B 2 F 2 Y 1 {r 2 BF 2 A' 2 + 2AF [2B 2 F 
+r 2 FA'B' - rB (A' (6F + rF') + 2rFA")) 
+A 2 (rFB 1 (4F + rF') + B (-4F 2 + r 2 F' 2 
-2rF{2F' + rF"))]} , (12) 

which involves first and second derivatives of the g a b met- 
ric coefficients F, B and A. By using the algebraic re- 



lations, we can write R g only in terms of the g a &-metric 
coefficients p and h and of the matter fields P and p. The 
final expression can be schematically written as 



R g = R g {p,p',p",h,h',P,P',P" 



(13) 



where we have used the equation of state to eliminate p 
and its derivatives. Note that in general relativity the 



Ricci curvature simply reads R g 



-T = p — 3P, i.e. no 



derivatives of the matter fields appear. This has profound 
implications. For example, if the function P(r) is contin- 
uous but not diffcrcntiable at the stellar surface, then P' 
would be discontinuous at the radius and P" would in- 
troduce an unacceptable Dirac delta contribution to the 
curvature. 

In order to see this point more clearly, we can use the 
field equations @ to eliminate the derivatives in Eq. f| 13[) . 
For simplicity, we shall assume an equation of state of 
the form P = A"/5q, where po is the rest-mass density 
and K and T are constants. The energy density can then 
be written explicitly as p{P) = [P/K] 1/r + P/(T - 1). 
Finally, by evaluating R g at the stellar surface, i.e. as 
r —> Rs and P — > 0, for any c^Owe get 



R g (P 




o < r < 3/2 
3/2<T<2 . (14) 

r > 2 



In the equation above, 7r = 0if0<r<3/2, whereas 



7r 



*(2-r) 
2ijA' 3 / r r 2 

-8k 2 



8(l-r)AT 1 / r 



3/2 < T < 2 

r = 2 
r > 2 



(15) 



with r] = R 3 s {Rs - 2M)/M 2 and M is the total mass 
defined by h[R s ) = [1-2A//As] _1 [SEl]. Therefore, for 
any T > 3/2 the scalar curvature diverges at the surface. 
The diverging terms originate from the derivatives of the 
matter fields in Eq. (|13p . 

A similar result, i.e. the divergence of the Ricci scalar 
at the surface of poly tropic matter configurations, was 
obtained in Palatini f(lZ) gravity and the consequences 
are discussed in detail in Refs. [l7l - fl9j ]. In that case, 
the divergence occurs only for 3/2 < T < 2. In EiBI 
gravity one has 7r ~ 1/k for T > 2, so that this singular 
solution does not appear in a small n expansion in which 
EiBI gravity resembles Palatini f(lZ,lZ a b). 

What we have established is that EiBI gravity does not 
admit any regular solution for polytropic spheres with 
r > 3/2, even for arbitrarily small values of k. At least 
two physical matter configurations are exactly described 
by a polytropic equation of state with T > 3/2: a degen- 
erate gas of nonrelativistic electrons and a monoatomic 
isentropic gas, both having V = 5/3. These perfectly 
reasonable systems, which can even be described within 
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Newtonian theory, have no description in EiBI theory. 
This renders the theory at best incomplete. 

It is worth stressing that our analysis relies only on the 
field equations and on the form of the equation of state 
close to the stellar surface. Any matter configuration 
whose behaviour resembles a polytrope (as an effective 
description) with adequate accuracy in the immediate 
vicinity of the surface will be singular, regardless of any 
complicated microphysics describing the stellar interior. 
There are many known examples of stars that satisfy this 
property. For instance the atmosphere of a white dwarf 
is well approximated by a polytrope with T = 10/3 (see 
Refs. 



e.g 
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At this stage one might claim that the polytropic, 
perfect-fluid description will break down at very small 
densities, and this might be a potential way out. 
Nonetheless, this would imply that no solutions for white 
dwarfs (and for many other systems) are allowed in EiBI 
gravity without precise knowledge of the microphysics of 
the matter near the surface. Even after introducing some 
microphysical description, there are no guaranties that 
the solution would be regular. Indeed abandoning the 
fluid approximation would just increase the differential 
order of the field equations in the matter sector, making 
the curvature even more sensitive to sharp variations in 
the matter fields. 

For the sake of the argument though, let us sup- 
pose that a polytropic equation of state provides a re- 
liable description close to the surface down to, say p ~ 
10 _n kg/ m 3 . Strong deviations from general relativity are 
expected when lBI S> Rf R = p — 3P at some radius 
very near the surface. For example, surface singularities 
would give rise to divergent tidal forces, which can be 
orders of magnitude larger than in Einstein's theory [l8[ . 
Let us then require that Rf lB1 < Rf R . This yields 



K > 4 x 10 24+2 ™m 5 kg-V 2 



(16) 



where we have assumed M ~ 1.4M , R$ ~ 10 _2 i?© and 
r = 10/3. For n = 10, the absence of strong near-surface 
curvature effects would imply n > 4 x 10 44 m 5 kg~ s~ 2 . 
This bound, though admittedly simplistic, is nonethe- 
less about 40 orders of magnitude larger than other cur- 



rent constraints [3j-l5|, [llj, [l4| . Additionally, Eq. (fT6"|) is 



a lower bound whereas all other constraints (including 
those one could derive by applying similar arguments to 
matter configurations well described by polytropcs with 
3/2 < T < 2) arc upper bounds that are totally incom- 
patible with Eq. (fTBj) . 

We emphasize that the surface singularities found 
above are not a prerogative of some polytropic fluid de- 
scription of matter near the surface. Polytropes just al- 
low for a convenient analytical treatment of the problem. 
The key issue is that higher-order derivatives of matter 
fields, which appear in the EiBI field equations as a re- 
sult of integrating out nondynamical degrees of freedom, 



make the geometry unacceptably sensitive to sharp vari- 
ations in the matter configuration. Surface singularities 
in polytropic spheres are just one possible manifestation 
of this sensitivity, but similar phenomena will be present 
in real world systems, where sharp density variations are 
common. This shortcoming, the rest of the pathologies 
related to the presence of nondynamical fields and a thor- 
ough discussion on the limitations of the polytropic fluid 
approximation can been found in Ref. [191 ] for Palatini 
f{1Z) gravity (see also Refs. 29|, 3(|). We will avoid re- 



peating this discussion here for EiBI gravity, as it would 
be nearly identical. 

In summary, we have shown that EiBI gravity is 
plagued by serious pathologies, whose root is the fact 
that the theory contains an auxiliary connection. The 
latter can be eliminated in order to obtain second-order 
dynamical equations for the metric only. However, the 
differential structure of these equations is profoundly dif- 
ferent from general relativity, as higher derivatives of the 
matter fields appear. This makes gravity noncommu- 
lative and spacetime geometry particularly sensitive to 
sharp changes in the matter configuration. In particu- 
lar, a discontinuity in the matter density or even only in 
its derivatives is enough to produce curvature singulari- 
ties, leading to unacceptable phenomenology. This fact 
has been missed in the recent literature on EiBI grav- 
ity but it has profound consequences for the viability of 
the theory, similarly to the case of Palatini f(lZ) grav- 



ities 
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and to theories where matter is coupled to 
the Ricci scalar [2oT ]. 

These problems appear to be a generic prerogative of 
gravitational theories which do not propagate any de- 
gree of freedom other than the massless spin- 2 field, but 
instead contain auxiliary fields that arc algebraically re- 
lated to the metric and to matter. Although these theo- 
ries are equivalent to general relativity in vacuum, once 
the auxiliary fields are eliminated using some algebraic 
relation, the dynamical gravitational equations contain 
higher derivatives of the matter fields. This is the case 
discussed here, where the auxiliary field is the connection 
or, equivalcntly, the q a b metric. 

In principle, this pathology could be alleviated by in- 
troducing higher order derivatives of the metric as well 
(this is for example the case in which matter is coupled to 
some nonlinear function of the Ricci scalar 311). How- 



ever, this usually produces different problems, like uni- 
tarity loss, instabilities, etc. 

Finally, in this work we followed the original proposal, 
assuming the independent connection to be symmetric, 
the matter action to be independent from it and that only 
IZ(ab) appears in the action. Relaxing these conditions 
would lead to the most general, metric-affine [25[ ver- 
sion of EiBI gravity, which presumably has a richer phe- 
nomenology. Understanding whether this more general 
theory may evade the pathologies discussed here would 
be an interesting extension of the present work. 
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